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Abstract 

We describe ellipticity domains for the isochoric elastic energy 


F !->• II dev„ log 17||^ 


VWf ^ 

(detF)i/" 


1 

4 


(detC)!/- 


for n = 2,3, where C = F'^F for F £ GL+(n). Here, devnlogH = log H - i tr(logl7) • 1 is 
the deviatoric part of the logarithmic strain tensor log U. For n = 2 we identify the maximal 
ellipticity domain, while for n = 3 we show that the energy is Legendre-Hadamard elliptic in 
the set 

£3(w^°,LR,U,^'^ :=|;7GPSym(3) | || devg log t/f < 11, 

which is similar to the von-Mlses-Huber-Hencky maximum distortion strain energy criterion. 
Our results complement the characterization of ellipticity domains for the quadratic Hencky 
energy VFjj (F) — jj, || devs log 17||^ + -^ [tr(log H)]^, U = VF'^F with /r > 0 and k > | /r, previ¬ 
ously obtained by Bruhns et al. 
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1 Introduction 

The quadratic Hencky energy 

(F) = /i II log C/f + ^ [tr(log U)]^ = II dev„ log t/f + f [log(det U)]^, 

where k > 0 denotes the bulk modulus and n, X are the Lame constants with /r > 0 and 3 A + 2 /r > 0, 
has recently been shown to have a fundamental geometric property which uniquely characterizes it among 
all hyperelastic formulations: it measures the (squared) geodesic distance of the deformation gradient F = 
€ GL~^(n) := {X G | detX > 0} to the special orthogonal group SO(n) [311155] . Here, U = 'JF’^F 
is the right stretch tensor and dev„logl7 = \ogU — ^ tr(logC/) • 1 is the deviatoric part of the Hencky 
strain tensor log C/, 1 denotes the identity tensor on R”^", ||Xp = {X,X) is the Frobenius tensor norm and 
tr(X) = {X, 1 ) is the trace oi X G R"^". The Hencky energy was originally introduced by Heinrich 
Hencky in 1928 [T5ll30| . see also Richter’s 1948 paper [ICTl eq. 7.3]. However, in an 1931 article in the Journal 
of Rheology m, Hencky also considered elastic energy potentials of the form 

Wig 3 i{F) = ^ II dev3 logHlI^ + h{detU). (1.1) 

Here, the volumetric part h : R_|_ := (0, oo) —^ R of the energy is a function to be determined by experiments. 
In an 1933 article [20], he suggested an even more general expression for describing the elastic behaviour of 
vulcanized rubber: 


1 ^i 933 (-F) = M ^ .f (^log Az A 3 )V 3 ) ^(detH) =/r/(dev 3 log[/) +/i(detC/), 

where / : Sym(3) —>■ R is an isotropic function in Valanis-Landel form |481 E5] . The ellipticity of W 1931 , 
provided that h is convex on R+, depends only on the ellipticity properties of the isochoric term || dev 3 log t/||^ 
measuring the purely distortional part of the deformation, which we investigate in this article. 

The necessity of finding an ellipticity domain for the isotropic invariant || dev„logC/|p of the logarithmic 
strain tensor log!/ (see |32l IMl [MlI42| ) arises from the observation that the isochoric part 

W'^°{F) :=M||dev„logC/||2 

of the quadratic Hencky energy is not rank-one convex even in SL(n) '■= {X G GL^(n) |detJf = 1} for 
n = 2 ,3 (see |37j). The understanding of loss of ellipticity is of fundamental importance in nonlinear elasticity 

[MlllsKlslIlIlEllEHKIilET]. 

It is easy to show that a given hyperelastic formulation is not rank-one convex. In general, it is also 
clear that there exists a neighborhood of the identity tensor 1 where the formulation is Legendre-Hadamard 
elliptic (LH-elliptic). What is difficult, however, is to precisely describe the maximal domain of ellipticity: 
although we are able to numerically determine the maximal ellipticity domain for || dev„ log Hp, in this paper 
we choose an intermediary way in that we analytically describe a large set in which the energy is LH-elliptic. 
The numerically obtained visualization in Fig. [2| indicates that this subset, expressed in terms of certain 
transformations of the principal stretches, is in fact the ellipse inscribed in the maximal ellipticity domain. 
For practical applications (like the coupling with elastoplasticity), knowing such a domain is mostly sufficient. 

The analysis in this paper is also motivated by the results established for n = 3 by Bruhns et al. IZIIHI 
(see also EU US] in order to compare the domains of ellipticity obtained in nonlinear elastostatics for a 
special material), who found an ellipticity domain for the quadratic Hencky strain energy W^: they showed 
that Wjj satisfies the Legendre-Hadamard condition for all principal stretches A^ with A^ G [0.21162..., ^fe\ = 
[0.21162..., 1.39561...], provided that the additional condition A > 0 holds. This result, however, is not 
applicable to the deviatoric quadratic Hencky energy || dev„ log C/p, which corresponds to the case A = 

-¥<o- 

It might also be worthwhile to find a scalar function of the isotropic invariant || dev„ log C/p such that the 
composition is elliptic over GL'''(ri). Indeed, in the two-dimensional case we have identified such functions 
EZldSl Eni|3S|, namely the so-called exponentiated Hencky energies 

Wl^{F) := F G GL+(2), (1.2) 

n, 


2 




where k is an additional dimensionless parameter. In m it is shown that these energies are polyconvex 
for fc > -j. In fact, (F) is polyconvex if and only if fc > |, see [2^, while || dev 2 logC/|p is not overall 
rank-one convex. For n = 3, however, such a function is not yet known. 

Knowles and Sternberg [22] have established a criterion for rank-one convexity (ellipticity) which is nec¬ 
essary and sufficient for n = 2 but only necessary for n = 3 (see also [2] [310 |T21 [231 |T| for alternative proofs). 
For n = 3, necessary and sufficient conditions for ellipticity were given by Simpson and Spector m, while 
for arbitrary dimension they were established for the first time by Silhavy [33] in terms of the copositivity of 
certain matrices. The necessary and sufficient conditions introduced by Dacorogna |9] were obtained by com¬ 
bining a result established by Silhavy |44| with one result on copositive matrices in dimension 3 by Hadeler 
m- In this paper, we use Dacorogna’s sufficient criterion |9| for arbitrary n, which can be applied more 
easily than Dacorogna’s necessary and sufficient conditions jS] Theorem 5], and for n = 2 is equivalent to the 
necessary and sufficient criterion previously shown by Knowles and Sternberg |22| . 

We will, however, not go into further detail on the general importance of rank-one convexity. For more 
information on this topic, we refer to the comprehensive books uniiisiiii] and to the papers lasiisisniin]. 

2 Preliminaries 

An energy W : GL“'' (n) —>■ R is called rank-one convex |3J page 352] on GL“'' (n) if it is convex on all closed 
line segments in GL''’(n) with end points differing by a matrix of rank one, i.e. if 

w{F + {i-e)^(^ri) <e wiF) + {i- e)w{F -t c g ??) ( 2 . 1 ) 

for all F G GL~^(n), 0 G [0,1] and all rj G R" with F + t^ G rj G GL“''(n) for all t G [0,1], where ^ G 77 
denotes the dyadic product. Using definition (12.IL in [29] it was shown for the first time that the mapping 
F I—>■ II dev„ log {7IP is not rank-one convex on all of GL~''(n). 

Since GL~''(n) is an open subset of R”’^”, an energy W : GL''’(n) —R of class is rank-one convex if 
and only if it is Legendre-Hadamard elliptic (LH-elliptic) at all points F G GL“'' (n): 

77 ,^ 077 ) > 0 for all 5 ,77 e R” . (2.2) 

Note carefully that, by this definition, rank-one convexity is strictly a global concept: a function on GL’''(n) 
is either rank-one convex or it isn’t. Legendre-Hadamard ellipticity, on the other hand, is also well-defined as 
a local property: a function W : GL^(n) —R is called LH-elliptic (or simply elliptic) on a set £ C GL“''( 7 t,) if 
( 12:211 holds for all F G £. In this case, £ is also called a domain of ellipticity or ellipticity domain for W. We 
also use the term maximal ellipticity domain to refer to the set of all points in which a function is LH-elliptic. 

Let us remark that whereas Bruhns, Xiao and Mayers ims] directly used definition ( 12 . 21 ) for finding an 
ellipticity domain of the quadratic Hencky energy, in this paper we do not calculate the second derivative 
DpW{F) of the energy W{F) = || dev„logC/|p. Instead, we will consider the representation of the isotropic 
energy W in terms of the principal stretches and utilize criteria applicable to this representation. 

Next we recall some of these useful results about LH-ellipticity as well as some properties of the deviatoric 
part of the strain tensor log U. 

2.1 Criteria for LH—ellipticity based on principal stretches 

In the three-dimensional case, our purpose is to identify an ellipticity domain, but not necessarily the maximal 
one, for the energy F 1 —>■ || devalogUlp. We therefore need a suitable sufficient criterion for LH-ellipticity. 
The following theorem was given by Dacorogna |2] Proposition 7] in the form of a criterion for rank-one 
convexity, i.e. for ellipticity on all of GL“''(n). It can easily be seen from his proof that the local form given 
here holds as well; note that the requirement that the set £ is open or the closure of an open set ensures that 
every F G £ can be written as the limit of a sequence {Fk)k C f of matrices with pairwise different singular 
values, which is utilized in Dacorogna’s proof. The criterion has previously been used by Gliige and Kalisch 
m in a similar way. 

Theorem 2.1. Let W : GL~'’(7t,) be an objective and isotropic function of class with the representation 
in terms of the singular values ofU via W (F) = g{Xi, X 2 , ■■■, A„), where g G C'^(R" , R) is symmetric. Further, 
let £ C GL’’" (n) be an open set or the closure of an open set. Define £ C R" as follows: 

(Ai,..., A„) G £ if and only if there exists F G £ such that Ai,..., A„ are the singular values of F . 
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Then W is Legendre-Hadamard elliptic at all F G £ if the following four sets of conditions hold: 


i) 



“TE—inequalities” 


for every 


* = 1, 2 ,n 


and all (Ai, A 2 ,A„) G f , 


(2.3) 


ii) for all i j , 


\ _ \ 
d\i ^3 dXj 

Xi Xj 

“BE—inequalities’ 

1 


> 0 


d'^g d^g d'^g 


n — 1 

1 

n — 1 


dXj dX] 


dXidXj 


Id^gd^g d^g 


dXj dX] 


dXjdX: 


da 

d\i 


da 

3A, 


Xi Xi 


da 

d\i 


da 

9A, 


Xi + A, 


> 0 


> 0 


For n = 2, the conditions are also necessary. 


for all {Xi, X 2 , ■■■, Xn) G £ with Xi ^ Xj , 


for all (Xi, X 2 , ■■■, Xn) G £ with Xi ^ Xj , 
for all (Ai, A 2 ,A„) e £ . 


(2.4) 


(2.5) 

( 2 . 6 ) 


Here M+ = (0, 00 ). The necessary and sufficient conditions of this theorem in the case n = 2 are the same 
as established by Knowles and Sternberg [23j|24], see also |43l page 318]. 

Dacorogna [9j page 6] also explains that due to the permutation symmetry of g, it is enough to establish 
only 4 inequalities: one TE-inequality (tension-extension inequality) for i = 1, one BE~inequality (Baker- 
Ericksen inequality) for i = 1 and j = 3 and two other inequalities from (12.51) . (12.61) for i = 1, j = 3. 
Note carefully that this remark is valid only when one considers the question whether a function is rank-one 
convex, i.e. LH-elliptic on all of GL“''(n); if a specific domain £ is considered, then the corresponding set £, 
which consists of all (Ai,..., A„) G R" which are singular values of some F G £, has to be invariant under 
permutations in order to reduce the number of inequalities. 

If, on the other hand, one wants to completely characterize the maximal ellipticity domain for an energy in 
spatial dimension n = 3, then the necessary and sufficient conditions of Dacorogna [5J Theorem 5] are better 
suited. In this set of conditions, one has to show 10 inequalities. However, due to some other symmetries 
and invariance properties, and since the BE-inequalities are always satisfied by || deVnlogf/|p, there remain 
5 inequalities in the necessary and sufficient conditions of Dacorogna which have to be checked in order to 
study the ellipticity of the energy || devalogt/p. We do not use this criterion in the analytic part of this 
article. 


2.2 Auxiliary remarks 

n 

The norm of the deviator in is given by || devn diag(^i,... ,^n)|P = ^ E ~ ■ Thus, for 

F S GL'''(n) with singular values Ai, A 2 ,..., A„, it follows that i,j=i,i<j 

II dev„logC/||^ = 5 (Ai,A 2 ,...,A„), (2.7) 

where the function g : M" —>■ [0, 00 ) is given by 

n 1 ” A 

5 (Ai,A 2 ,...,A„) = - ^ (logA* - logA jE = - ^ (2-8) 

Tl TL /\j 

Note that the function g is invariant under scaling: 

g{a Xi,a X 2 , ■■■,a Xn) = g{Xi, X 2 , ■■■, Xn) for all a > 0. (2.9) 

Hence, for the function g corresponding to our energy E || devji log [/|p, the inequalities in Dacorogna’s 
criterion are also invariant under scaling, see m for further details. Therefore, for an arbitrary scaling factor 
a > 0, the function g satisfies the required inequalities from Dacorogna’s criterion (^Theorem 12.II) in a point 
(Ai, A 2 ,..., A„) = (a Ai, a A 2 ,..., a A„) if and only if it satisfies them in the point (Ai, A 2 ,..., A„). 


4 












Since || dev„ logJ7|p linearizes to ||dev„e||^, where e denotes the linearized strain tensor, it is obvious 
that the maximal ellipticity domain of || devji log f/|p contains a neighborhood of 1. Moreover, the above 
considerations show that this domain is an (unbounded) cone containing 1. In the following we will exploit 
this insight. 


3 The two-dimensional case 


Using the ellipticity conditions by Knowles and Sternberg, i.e. Theorem 12.11 for n = 2, we obtain: 
Proposition 3.1. The maximal ellipticity domain of the energy T i—>■ || dev 2 logC/|p, F € GL“'"(2) is 

:=|c/ePSym( 2 )|||dev 2 logt/f < i|. (3.1) 

Proof. We will prove this result using the necessary and sufficient conditions given by Theorem 12.II for n = 2 
together with the identity (12.711 . To this aim, we need to compute 


dg 1 1 ^1 

— = — log — , 

uXi Ai A2 


dg 1 1 -^1 

WT = “W 
aA 2 A2 A2 


(3.2) 


^ = i 

dXj Xl 


- log ^ + 1 


d^g 


1 


^2 L 


log ^ + 1 
A 2 


d^g 

dXidX2 


Xl X 2 


and verify that inequalities (I2.3I) - (I2.6I) hold if and only if (Ai,A 2 ) G £ 2 , where the set £2 of singular values 
corresponding to the domain £ = £2 (lUg®°, LH, U, i) is given by 


log’ (i < 1 

A 2 


'£’2 = (Ai, A 2 ) G K 

The TE-inequalities of Theorem 12.11 are equivalent to 

- log ^ + 1 > 0, log ^ + 1 > 0, (3.3) 

A 2 A 2 

while the BE-inequalities are satisfied everywhere for convex functions of log U m and thus in particular by 
II dev 2 logC/|p. The inequalities (12.41) and (12.51) are equivalent to 


- log ^ + 1 
A 2 


log ^ + 1 
A 2 


- log ^ + 1 
A2 


log 


Ai 

A 2 


1 


- 1 + log ^ > 0 if Ai ^ A 2 , 

Ai — A2 A2 

Ai — A 2 -^1 ^ n 


(3.4) 


(3.5) 


Since all these inequalities are symmetric in Ai and A 2 (and thus the ellipticity domain is invariant w.r.t. 
the transformations Ai 1 —?► A 2 , A 2 1 —>■ Ai) we may assume that Ai > A 2 , i.e. that t '■= ^ > 1. Geometrically 
speaking, considering this substitution means that it is necessary and sufficient to prove that the inequalities 
(ixa . (HI and (13.51) are satisfied along all lines Ai = t X 2 , t > 1. 

Thus, the inequalities (13.31) . (13.41) and (13.51) are satisfied if and only if the following inequalities hold: 


-logt + 1 > 0 , 


logt + 1 > 0 , 


\/[-logt + 1 ] [logf + 1 ] - (1 - logt) + -—- logt > 0 


t - 1 


\/[-logt + 1 ] [logt + 1 ] + (1 + logt) + ——- logt > 0 . 


t 1 


(3.6) 

(3.7) 

(3.8) 


Since t > 1 and y/[— log t + 1] [log t + 1] > (1 — logt) for all t > 1 , we find that the inequalities (13.71) 
and (|3.8I) are redundant in the set of inequalities describing the domain of ellipticity. In conclusion, the 
independent inequalities describing the ellipticity domain are 

— logt + l>0, logt + l>0, (3.9) 

which can equivalently be expressed as 1 < log^ t = log^ ^. Therefore we deduce that the ellipticity conditions 

are satisfied if and only if (Ai,A 2 ) G £ 2 , i-e. if and only if || dev 2 logt/|p = 5 log^ < 1, and the proof is 
complete. ■ 
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4 The three-dimensional case 

For n = 3, we consider the substitution 


Then 



^ ^ ^ p-(a+&) 

As ’ Ai 


(4.1) 


Ai _ A2 _ As 

^ “ T “ ^ 


(4.2) 


which means that (Ai,A 2 ,As) belongs to the line which passes through ( 0 , 0 , 0 ) and an arbitrary point 
(e“,l,e“^) in the plane A 2 = 1. According to the invariance properties of the energy and of the condi¬ 
tions for ellipticity given in our preliminaries, it is enough to study the resulting inequalities only in the plane 
A 2 = 1. 

Numerical calculations indicate that the three-dimensional maximal domain of ellipticity is that for which 
(a, b) in (14.11) belongs to the two-dimensional domain described in Fig. [TJ However, since it is difficult to 
characterize the maximal ellipticity domain, we consider a significant large subdomain of it (see Fig. 




Figure 1: The maximal ellipticity domain 
in terms of (a, b) obtained numerically after 
implementation of the necessary and suffi¬ 
cient criterion of Dacorogna [HI Theorem 5], 
which we do not use in the analytic part of 
this article. 


Figure 2: We consider the full ellipse which 
is contained in the maximal ellipticity do¬ 
main and which touches its boundary. 


Proposition 4.1. The maximal ellipticity domain of the energy F 1 —>■ || deva log t/|p, F £ GL~'’(3) contains 
the ellipticity domain 

fa (^Wr,LH,t/,0 := |[/S PSym(3) | ||deva log t/f < || . (4.3) 

Proof. According to 12. II we have to show that inequalities (I2.3I) - (I2.6I) hold for all (Ai, A 2 , Aa) £ S 3 , where the 
set S 3 of singular values corresponding to the domain S 3 (11T®°, LH, U, |) is given by 

fs = |(Ai,A 2 ,Aa) £ M+ log^ ^-h log^ ^-f log^ ^ < 2 

A 2 A 3 A 3 

Let us first observe that, under the substitution dm, 

; -^2 ; -^ 3 ) ^ ^3 ^ + 6^ + a 5 < 1 , 
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which means that (a, b) belongs to the full ellipse + ab = 1. This domain is invariant under the three 

transformations 


a = x + y, b = -y] 


a = y, h = x] 


a = —y, b = —X (4.4) 


in the sense that if {x, y) belongs to the full ellipse w'^ +z'^+w z = 1 , then (a, b) belongs also to this domain and 
vice versa. This invariance property is equivalent to the (obvious) invariance of the set £3 under permutations 
of Ai, A 2 and A 3 . From this symmetry, it follows that we only need to consider the inequalities in Dacorogna’s 
criterion for i = 1 and j = 3. We compute 


^ = J_ A ^ 

5Ai 3Ai A 2 
d^g _ 2 / A 3 

d\l 3 A? (, °®Ai 

The TE-inequality for i 


A 3 

Ai 

Ai 


- log 

- log 
= 1 is equivalent to 


0 X 3 

d'^g 


dXl 


2 

3 A 3 
2 


1 A 2 1 A 3 

log ^ ^ 


-log 


As 


Tog 


A 2 

As 


2 “h ct — b ^ 0, 


(4.5) 

d^g _^ 

8 X 38 X 1 3 A 1 A 3 ■ 


(4.6) 


while the BE-inequality for i = 1, j = 3 reads 

A 3 3 ^ — Ai gj;- _ 1 —2 (g + 2 6) — 2 (2 g + 5) _ 1 —2 (a + b) ^ 
A 3 - Ai ~ - 1 ~ W 6-1“+'’) - 1 “ 

and is always satisfied. Moreover, we compute 


and 


1 j 8 ^g 8 ‘^g 

2 Y 5Ai 8 X 1 


8 '^g 


8X38X1 

2 

3 A 3 Ai 
2 

3 A 3 Ai 


T 

.2 

T 

.2 


dg I 

0X3 ' dXi 
As + Ai 


a/( 2 + a + 2 6) (2 - 2 a - 6) + 1 + 


v/(2 + a + 26)(2-2a-6) + l + 


—Ai (a + 2 6) + A3 (2 a + 6)' 

A 3 + Ai 

-e“+'' (a + 26) + ( 2 o + 5)] 

1 + e“+'' 


1 / 8 ‘^g 8 '^g 

2 Y 


02^ ^ ^ 

O g , d \3 d\i 

8X38X1 A3 — Ai 


2 

3 A 3 Ai 
2 

3 A 3 Ai 


^ 3/{2 + a + 2 b){ 2 - 2 a-b)- 1 + 
i a/( 2 + a + 2 6) (2 - 2 a - 6) - 1 + 


—Ai (u T 2 6 ) — A 3 (2 Q b^ 

A 3 — Ai 

-e°+fa (a+ 2 6) - { 2 a + b)- 
1 — e“+^ 


We therefore need to show that for all a, & € R with + b^ + ab < 1, the following inequalities, corresponding 
to conditions (USD, (EH) and (12.61) for i = l,j = 3, hold: 


— 1^(2 + a + 2 6) (2 — 2 a — 6) + 1 + 
^ 3/(2 + a + 2b) {2 - 2a - b) - 1 + 


2 + a-5>0>0, 
-6“+''(a+ 2 6) + (2a + 6) ^ 

1 + e“+^ - " ’ 

-6“+''(a+ 2 6) - (2a + 6) ^ 

1 - e“+& - " 


if a + b ^ 0 . 


(4.7) 

(4.8) 

(4.9) 


Again, explicitly writing out the required inequalities for all i,j S {1,2,3} with i ^ j would simply yield 
inequalities which can be transformed into (I4.7I) - (I4.9I) via the transformations (14.41) . Using the further sub¬ 
stitution 


a = 


(a/s cos 0 -I- sin0) = p 



b = -^(—a/ScosS -|- sin0) = —p 
V6 
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Figure 3: Graphical representation of the functions /i,/ 2,/3 defined in (14.121) . 


we find 


and 


2(a^ + 5^ + ab) = 2 ^(3 cos^ 9 + 3 sin^ 9) = 


. , ‘^P ■ a 

a + b = —j= Sinn, 
V6 
P 


. — b = \f2p COS 0 , 


a + 26 = v^cosfl + 3 sin0) = ^PXp sin (j) — —^ 

TT' 

6. 


2 a + 6 = (a/s cos 0 + 3 sin0) = V^p sin ^0 + . 


(4.10) 


(4.11) 


These substitutions imply that the point (a, b) lies on the ellipse 2(a^ + b^ + ab) = p^. Note as well that for 
all p G R the corresponding ellipse is invariant under the transformations (14.41) and that (a, b) lies inside the 
full ellipse + 6^ + a 6 = 1 if and only if p < V^. In the following we prove that the condition p < a/^ is 
sufficient for ellipticity. 

In terms of p and 9, the required ineqalities (I4.7I) - (I4.9I) can now be written as 


flip, 9) :=2 + V2p COS 6* > 0 


f2ip,9):=^^ 2 + V2psm(9-^^ 2 - v^p sin (^0 + 0 


^ ^ ^ V2psin(0-f)+y2psin(0+§) ^ ^ 


-I , sin 9 

1 + 


fsip, 2 + sin 2 - \/2p sin ^9 + 0 

_ ^ ^ -e7^"‘"^V2p sin (0 - - V^p sin (0 + g) ^ ^ 


I — qVs 


sinfl 


V(p,0)G [0,V^] X [0,2 7r), 


(4.12) 


V(p,0)G [0,^]x [0,2 7r), 


V(p,0)G [0,v^]x{(0,2^)\{7r}}. 


We observe that the functions /i, / 2 , /a are periodic in 0 with period 2 tt. Moreover, from 
sin(2 TT — 0) = — sin 0 , cos(2 tt — 0) = cos 0, 
sin ^2 TT — 0 — —^ = — sin + 0^ , sin ^2 tt — 0 + —^ = — sin ^ (4-13) 

it follows that 

flip, 0) = flip, 2 71 - 9), /2(p,0) =/2(p,2 7r-0), fsip, 9) = fsip, 2 77 - 9). (4.14) 


It therefore suffices to show that the functions fi, /2 and fs are non-negative in the domains [ 0 , a/2] x [0 ,7r], 
[ 0 , a/2] X [ 0 , 7 r] and [0,a/2] x ( 0 , 7 r), respectively. 

These functions are indeed non-negative on these domains, as the graphs in Fig. clearly show; the 
reader can find a complete analytical proof in Appendix lA.ll This last assertion completes the proof of our 
proposition. ■ 
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Figure 4: The obtained ellipticity domain 
of the energy || deva logf/|p in the principal 
stretches Ai, A 2 , A 3 . 


Figure 5: The intersection of the ellipticity 
domain with the plane Ai + A 2 + A 3 =3. 


5 Concluding remarks 


In order to visualize the established domain of ellipticity for u = 3, we go back to the initial substitution and 
find that the unbounded ellipticity domain given by Proposition 14. II is the set enclosed by the cone presented 
in Figures 0] and [S] which is completely defined by 

Ai X 2 =u, A 3 =ue*(^“^“'), 0 e [0,2^), u € [0, 00), p e [0, v^]. 


It is clear that the ellipticity of the energy T 1 —>■ p || deva log C/|p on GL’^(3) for /i > 0 in the domain 


£3 LH, [/, - j = |c/ e PSym(3) | || deva log Gf < 

implies the ellipticity in this domain for the exponentiated Hencky energy 


^ ^ gl|dev 3 log (7|p gfc [tr(log U)f 

2 k 


p, K > 0 , k > -, 


(5.1) 


(5.2) 


since on the one hand t 1-^ fie* is convex and monotone increasing, and therefore the composition with this 
mapping preserves ellipticity, and on the other hand the function F 1 —>• gk (logdet f)'^ jg j-ank-one convex on 
GL“'"(3) for /c > i (see [37] for more details). 

However, numerical tests suggest that the ellipticity domain of the exponentiated Hencky energy is far 
bigger than all ellipticity domains which are known for various energies of quadratic Hencky energy type, see 
also |37l[36]. This remark might be useful in the study of large deformations which do not belong to the 
known ellipticity domains of the energies of the quadratic Hencky energy type. 

Note that the ellipticity domain £3 (W))™, LH, G, |) conforms exactly to the von-Mises-Huber-Hencky 
criterion, also known as the maximum distortion strain energy criterion in elasto-plasticity. Based on the 
results of the present paper, it is clear that the quadratic Hencky energy coupled to multiplicative plasticity 
will never lose LH-ellipticity in elastic unloading. This claim has been detailed in |36l l35l [34] . We do not 
know of any other elastic energy in which the ellipticity domain and the elastic domain coincide in this way. 
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A Appendix 


A.l The positivity of the functions fi, f 2 , fs in Proposition 14.Il 

It remains to show that the functions /i, / 2 , /s from the proof of Proposition 14.11 are non-negative on the domains [0, y/2] x [0, tt], 
[0, \/2] X [0,7r] and [0, \/2] x (0, tt), respectively. The condition fi{p,0) > 0 simply reads 

\/2 > —p cos 6 V 0 E [0, tt] ('^•3) 

and is obviously satisfied for all p E [0, \/2]. In order to prove that fsip, 0) > 0 for all p E [0, \/2] and 9 E (0, tt), we will prove 
more: we will show that the inequality 


sin(0- 2.) -V2psin(6»+f) 

-1 +- - —i^ > 0 

—£ sin 9 

1 - e V6 


(A.4) 


holds for all p E [0, y/2\ and 6 E (0, tt). Note that for {p, 9) E [0, ^/2\ x (0, tt) we have e since sin0 > 0 for 9 E (0, tt). 

Therefore, instead of proving the inequality IIA.4II . it is enough to show that 


1 — e V6 ^ -I- \/2pe ^ sin ^ + V2p sin {9 -h —^ > 0 


2 p • 2 p 

—C. cm H ,— —C. ^ , . 

Sin [ 9 — , , . _,_, - , 

" V 6^ 

for all p S [0, \/2] and 9 S (0, tt) . To this aim, let us remark that for 6 £ [ ^, tt) , we have 

sin 6 /— sin (? . f r, /iT -in \ 

1 —ev^ +\/2pev^ sm (^6 — — j + \/2p sm + — j 

>1 _ e ^ + \/2p sin ^6- ^ + \/2p sin - ^ =1 — + \/6p sin 9 > I 


Here we have used that 


sin (9 + —'i + sin (9 — —= 2 cos — sin 6 = \/3 sin 6 

V 6/ V 6/ 6 


(A.5) 


(A.6) 


(A.7) 


and that the function f 1 —> 1 — e ^ + \/6C non-negative on [0, \/2]. This shows that f3{p,9) > 0 for all p £ [0, \/2] and 
9 e [f ,vr). 

Let us now consider the inequality llA.4b for all p E [0, \/2] ^'iid 0 E (O, ^). We note that 

^/3 


, 7r\ \/3 1 

sin [9 -= -sin 9 -cos 9, 

6/2 2 


/ 7r\ y/S 1 

sin [9 -\ -) = -sin 9 -\— cos 9, 

V 6/ 2 2 


(A.8) 


11 





and that 9 £ (O, ■!■), cos9 > 0. We introduce the substitution C — P sinO, r] = p cos9, and our new aim is to prove that 




2 2 


+ -e^ +1>0 


(A.9) 


for all 77 £ [0, \/2] and £ [0, v^]- Again, this last inequality is stronger than necessary, since we are only interested in the case 
< 2, i.e. in a subdomain of [0, \/2] x [0, ^/2]. Inequality llA.911 is satisfied since 




ac 


Vs Ve 


(A.IO) 


> min 

C,i)e[o.V2] 


|eVl« ^ = ^J^>0 V(C,r;) S [0, ^2] X [0, ^2] 


implies that 

= 0 V(C, ^ 7 ) £ [0, \/2] X [0, v^]. (A.11) 

Thus the function f3{p,0) is also non-negative for all p £ [0, v^] and 9 £ (O, ^). Combining this with the earlier result for 
9 £ [^, tt), we find fsip, 9) > 0 for all (p, 9) £ [0, \/2] X (0, tt). 


In a similar way we remark that 


is equivalent to 


for all 9 £ [ 0 , 7 r]. Using that 


we find 


—e ^/2p sin (9 — J-) + V^P sin (^ + 5) 

1-h- - -^- — > 0 

—i sin 9 

l + e^ 


1 + V2p sin (e + J) + [l - sin (d - ^)] > 0 

sin (9 H-— sin (p -^ = 2 sin — cos 9, 

\ 6 '' ^ 6 6 


Ve6|0,-J. 


sin ((?+ — ) > sin ~ 


(A.12) 

(A.13) 


Thus, for all 9 £ [O, and all p £ [0, V^], we have 
1 + V2p sin (j) + —) + e ^ |^1 — V2p 


sin (6 — — )| + + \/2 p sin ( 6 - 


1 _ e ‘ 


, -S^sinS VG . „ 

= 1 + e H- p sin 9 

2 

^ sinS Vs . „ 

> 1 + e H- p sin 9 

2 


1 — eVe " 


sin 0 

1 - e 


- p COS 9 

2 


> 0 , 


sin 9 

1 - 


is positive on [ 0 , \/ 2 ], where 


Z P Q Z _ Z i, 

since cos 9 > 0 and 1 — e v% < 0 on [O, and since the function ^ 1 + e ( 1 — 

the substitution C = P sin0 was considered. The above inequality shows that / 2 (p, 9) > 0 for all p S [0, \/2] and 9 £ [O, . 

Next we prove that f2{p,9) > 0 for 6 S (^, 7 r]. We find 


1 + e 


h{pV) ^l + e'/e y[2 + \/2psin(6 »-[2 - \/2p sin ^0 - 

+ + -e^=“®V 2 p sinf 0 --)+V 2 p sinf 0 +-), 


1 sin 9 

= - ( l + ev^ 


4 - 2p2 sin^ [9 -)+ l + ev^ 


sine 


(A. 14) 


for all p S [0, \/2] and 9 £ ( ^ 5 tt] 5 since 


- “"®V 2 p sin {e-^)+ V2p sin (0 + 

in( 0 +n <sin(e-^) V 0 e(J,.]. 


(A.15) 


Combining iTOli with the fact that cos 0 < 0 for all 0 G (5 ’ introduce a new substitution C — P sin0, zo = —p cos and 

our new goal is to prove that 


/i(C,ro) ( 1 -he^/6 


4-2 ( ^C- iro) +(l + eVS 


Vs, 1 


— e^V2 ( — C + - 137 I + V2 ( — f 


Vs, 1 


= 2 ( + 1 ) e(C,ro) > 0 


(A.16) 
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for 3<11 ^5 uj G [0, v^], where 


I - ^ 

/ 3C^ zu^ e — 1 

s(C, ro) := W---h V3(u7 -— + 4 — \/2ro + 2 — \/6 C —- 

V 2 2 eV5+l 


(A.17) 


Again, this is more than is needed, since for the non-negativity of /2 for p G [0, ^/2] and 9 G (^if is enough to prove that 
vj) > 0 only for all ru E [0, \/2] which belong to the smaller domain < 2. We observe that 


—( 
d'uj 


Hence, we deduce 


1 

f 


(i 

( 

73C - ro H 

72 

U 

^-3C2 + 273Cro - ro2 + 8 j C,“6[0,v^ 

72 

VV 

-h 2\/3C ^ ~ + 8 y J 


^3 - 2 


< 0 VC,ro S [0, Vi], 


h(Ct^) > ^(C,v^) = \l + \/6C + 3 — 


Ve ^ 


+ 1 


(A.18) 


> max /i(C, V2) = f \/2-2-^3tanh (All Ri 0.0573242 > 0 VC,roe[0,V2], (A.19) 

Ce[o,\/2] \ \v3// 


which means that h is non-negative on [0, \/ 2 ] X [0, \/2]. This last conclusion shows that / 2 (p, 0) > 0 for all (p, 8 ) S [0, \/2] x(f ,7r]. 
Therefore, the function /2 is non-negative on [0,^2] x [0, tt]. 
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